The odd-even effect of the mehing temperature of polymer film on finite lattice 
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In order to find a quantitative understanding on the odd-even efi'ect of melting finite polymers, we 
compute the melting temperature and entropy growth rate of confined dimer film on finite rectangle 
and torus. The theoretical melting temperature demonstrates similar behavior as the experimental 
observations. For an infinitely long dimer film belt with finite width, the melting temperature 
strongly depends on the odd-evenness of the width. The entropy for an even number of width is 
always larger than the entropy for an odd number of width. When the length of the rectangle 
goes to infinity, the speed of entropy growth shows a finear dependence on the width. This linear 
relationship holds both for rectangle and torus. Fusing two small rectangles with odd number of 
length into one big rectangle gains more entropy than fusing two small rectangles with even number 
of length. Fusing two small toruses with even number of length into one big torus reduces entropy 
instead of increasing it. While fusing two small toruses with odd number of length would increase the 
entropy. The entropy difference between covering torus and covering rectangle decays to zero when 
the lattice size becomes infinite. The correlation function between two topologically distinguishable 
loops on torus also demonstrate odd-even effect. 

PACS numbers: 1.2.3.4.5.6 



I. INTRODUCTION 



Chemical studies found the melting point of fatty acids 
with odd number of chain length are below those with 
even number of chain length[l|. This odd-even effect 
was clearly shown by the melting temperature of the n- 
alkanes C„i/2n-i-2[l]- When the length of the n-alkanes 
increases by one unit from an even number to its next 
neighboring odd number, the melting temperature in- 
creases by a very small value compared with case for odd 
number growing to even number QQ. The odd-even ef- 
fect is diminished when the chain becomes longer. The 
miscibility of two different compounds also exhibit odd- 
even effect. Two different compounds dissolve in each 
other only when the length of the linker between them is 
odd [3. The odd-even effect is qualitatively understood 
by considering how to pack large molecules with differ- 
ent structure. A quantitative understanding of odd-even 
effect is still an open question, especially when it comes 
to anomalous odd-even effectQ. 

Dimer is the shortest polymer. Classical dimer model 
counts the number of different ways to cover a graph by 
pairing up two neighboring vertex[5j. The exact number 
of dimer covering on square lattice can be calculated by 
Kasteleyn's methodQ. The chemical polymer usually is 
longer than a dimer. However there is a theoretical corre- 
spondence between dimer configuration and fluctuating 
string configuration on two dimensional lattice [7|. Soft 
long polymer may bend into different shapes due to fluc- 
tuating environment. So we can get a basic theoretical 
understanding on melting long polymer by studying the 
melting of dimer film. 

The article is organized as following: In the section II, 
we computed the melting temperature of dimer film on 
a constant area. The shape of the area varies from even 



to odd or vice versa. In the section III, we take an in- 
finitely long dimer film with finite width, and study how 
the entropy growth rate behaves when we increase the 
width. In section IV, we study how two small rectangles 
fuse into one big rectangle. The same procedure for com- 
bining two small torus is computed to compare with the 
case of rectangle. In section V, the correlation function of 
two loops on torus is computed to check its dependence 
on odd-even effect. In section VI, we proposed the gen- 
eral framework for calculating the melting temperature 
of long polymers. The last section is a summary. 



II. THE MELTING TEMPERATURE OF DIMER 
FILM ON FINITE RECTANGULAR LATTICE 
WITH CONSTANT AREA 

The melting temperature of polymer is the critical tem- 
perature when a crystalline of polymers transforms into 
a solid amorphous phase. Dimer is the simplest poly- 
mer. We assume the crystalline phase of dimer film is 
one frozen pattern of the dimers. Each dimer is locally 
fixed in one direction. They can not rotate around. In 
the solid amorphous phase, the weak external bond con- 
necting neighboring dimers is broken. The dimers are 
free to rotate and reorient itself. But we can not make 
sure it must be in certain direction. So all different con- 
figurations are possible. 

We assume this phase transition is a first order tran- 
sition, then the Gibbs free energy, AG = A[/ — T/S.S, is 
zero at the phase transition point. We need to calculate 
how much thermal energy the film has absorbed to break 
the external bonds and the entropy of the two phases. 

Suppose the dimers are frozen in the beginning. When 
the dimers absorbed the energy of external heat source, 
it becomes active and oscillating around local equilib- 
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rium position. But the internal bond connecting two 
monomers is strong enough to survive during the mehing 
process. At very high temperature, dimer may also de- 
compose into monomers, but that temperature is much 
higher than the melting temperature. We denote the 
energy quanta for breaking a frozen external bond as 
Si. The total thermal energy for melting the film is 
AU = J2i ^i- We denote the partition function of the 
frozen phase as Zi, and the partition function of the 
melted phase as Z2. Then the critical temperature can 
be computed by 



AU 

'as 



ks log[Z2] - ks log[Zi] ■ 



(1) 



We consider a finite rectangle of square lattice with m 
rows and n columns. The total number of lattice sites is 
mxn. The total number of bonds is {2mn — n — m). The 
total number of internal bonds within a dimer is mn/2. 
The total number of external bonds is |mn — n — m. The 
thermal energy would shake the external bonds, when 
all of the external bonds become weak enough to reach 
the threshold of collapsing point, the dime film melts. 
The dimers randomly rotate to form a new pattern. The 
thermal energy absorbed by the dimer film at the critical 
temperature is AU = {^mn — n — m)eo. €0 is the unit 
energy of each external bond. We take it as eo = 10^^^ 
for the convenience of computation. 

The entropy of the frozen phase is zero since the num- 
ber of configuration is just one. In the melting phase we 
have to count all possible dimer coverings on the rectan- 
gle lattice. Kasteleyn had developed a method to count 
all possible ways for covering the whole lattice by dimers 
without isolated sites First one needs to construct 
Kasteleyn's adjacent matrix K by assigning each bond 
of square lattice a number -1-1 or —1. The product of 
these numbers around an arbitrary rectangle must be —1. 
The total number of all p ossible configurations of dimer 
covering is Z = \J det K [6|. Kasteleyn's method was 
equivalently mapped into quantum field theory of free 
fermion model [8| . A physical picture of dimer covering 
is to place a fermion at each lattice site. Two neighbor- 
ing fermions can form a pair in different direction. There 
is a hard-core repulsive interaction between two dimers 
since they do not overlap each other. 




FIG. 1: The dimer fihn is almost infinitely long 
varies from 2 to 16. 



The width 



We computed the meting temperature of dimer film on 
a large rectangle. The total number of covering lattice 
sites is mn = Uj^-^i x 11 x 13 = 6.16287 x lO^^. We 
keep this number invariant but varies the ratio of height 



to width. We let n running from 2 to 16. The length of 
this rectangle at n=2 is (6.16287/2) x 10^^, it is approxi- 
mately ten times longer than the length at n=16. When 
the width grows from 2 to 16, the dimer film grows from 
a long belt to a short belt with increased width. It can 
be viewed as folding an almost infinitely long polymer. 
The melting temperature increases as the width grows 
larger(Fig. [5]). When the width grows larger than 16, 
the melting temperature approaches to a constant num- 
ber. Then it reaches the thermal dynamic limit, the finite 
size effect disappeared. The increase of melting temper- 
ature at each step depends on the odd-evenness of the 
width. For example, when the width grows from 2 to 3, 
the critical temperature increased by 0.175617 percent. 
While it only increased 0.0561494 percent from 3 to 4. 
When we observe the whole curve of ATc, it shows the 
slop of ATc from an even number to an odd number is 
always larger than its sequent neighboring case from an 
odd number to an even number. This theoretical behav- 
ior of melting temperature is quite similar to the melting 
temperature of n-alkanes C„-ff2n-i-2[2|3- 
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FIG. 2: The melting temperature of dimer film at different 
width from 2 to 16. 
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FIG. 3: The increase of melting temperature at each step 
where the width of the dimer film grows bigger by one unit. 
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FIG. 4: The entropy of dimer coverings on the infinitely long 
rectangle. The total lattice sites of the rectangle is mn^llj^ii 
xll X 13=6.16287 x 10^**. The width of rectangle is repre- 
sented by m. 

At the melting point, the dimers are not strongly con- 
fined in a local position. Some dimers maybe jump out 
of its home to another hole. This provides some proba- 
bility to reshape the soft film. We can assume the long 
rectangle can vary its own shape following the second 
law of thermodynamics. When the rectangle reshapes it- 
self, it would exert some force to its surroundings. We 
call this force entropy force. As the shortest length 
n = (6.16287/16) x 10^^ is almost infinitely larger than 
the longest width m = 16, we view the length as infinity. 
We define the entropy force induced by increasing m as 



T 



AS'(to, n) 
Am ' 



(2) 



T is temperature. Here we take it as a positive constant 
for simplicity. Numerical computation shows the entropy 
at even number is larger than that at odd number(Fig. 
|4]). If Am = 3 — 2, one can see AS'(m,n) < by direct 
observation of Fig. HI So entropy force /(n, m(2 — >• 3)) 
is negative. For Am = 4 — 3, AS{m, n) > 0, the en- 
tropy force f{n, m(3 — > 4)) now becomes positive. This 
positive entropy force is much larger than the negative 
entropy force. If we enlarge the step size, Am = 4 — 2 = 
6 — 4 = • • • , the entropy force is constantly positive. Thus 
the direction of entropy force depends on the unit scale. 
When n grows larger, the entropy difference between even 
number and odd number vanished. So the entropy may 
experience periodic decreasing in finite sized system be- 
fore it reaches the maximal point. 



III. THE MELTING TEMPERATURE OF 
DIMER FILM ON A GROWING RECTANGLE 

We focus on the entropy increasing when the number 
of columns increases one by one. The number of rows is 
kept at a small value. Here the entropy of finite system 
is defined as the same formulation of thermodynamic en- 
tropy, 5* = kB^og[Z], Boltzmann constant ks is set to 




FIG. 5: The width of the dimer film is finite and increasing 
step by step. The length of the dimer film would be infinitely 
long. 

1 for convenience. For fixed number of rows, we com- 
pute how much entropy increased when the number of 
columns is increased by one unit. 



AS" = log [Z(m, n)] - log [Z(m, n - 1)]. 



(3) 



Fig. |6] shows the entropy increase on a rectangle and a 
torus when the number of columns increases from 2 to 
20. The number of rows is m = 10. When the length 
grows from odd number to even number, the step size of 
entropy increase is much larger than the case of grow- 
ing from even to odd(Fig. As the number of length 
grows, the entropy increase decays at even number, while 
the amplitude of AS' increases at odd number. This phe- 
nomena holds both for rectangle and torus. Different 
boundary condition does not smear out the odd-even ef- 
fect. 

When the length becomes longer, both the two 
branches of AS for even number and odd number con- 
verge to the same limit value. It takes more steps 
for AS" to converge on a torus than that on a rectan- 
gle. The limit value AS{n — > oo, m = 10) on rectangle is 
ASrecin oo, m = 10) = 2.778441, here I have made a 
cut off at 10^^. The limit value at large n on a torus 
is a little bit higher, AStor{n ^ oo,m = 10) = 2.969359. 
The difference between the two limit values of rectangle 
and torus preserves until m oo. Fig. [7] showed the 
limit values of AStorin oo) and ASrecin — oo) when 
m grows from 2 to 36. Both the two limit values fall on a 
straight line. For large numbers of m, m > 100, the limit 
values can be fitted by empirical linear equations. 



ASrecin 0O,m) = 0.29 m + Crec(™), 

AStorin ^ oo,m) = 0.29 m + Ctor("^), 
m ^ 2k, fc = 0, 1,2,3, 



(4) 



Without losing the key features, here the accuracy of 
those numerical values in equations are kept to the order 
of 10~^ for brevity. This empirical linear equation implies 
a constant. 
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An 



0.29. 



(5) 



This constant maybe has some hidden relationship with 
the entropy per site. However k deviates far from the 
exact number of the entropy per site in Ref. 
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ASrec(n — OO) 


AStor{n -> oo) 


^Stor ^Srec 


10 


2.778441055719985 


2.969359257199864 


0.190918 


20 


5.688333696422393 


5.857553050376143 


0.169219 


30 


8.601926338829257 


8.764325608871156 


0.162399 


40 


11.51650404237699 


11.67554509633454 


0.159041 


50 


14.43148654142293 


14.58852686413148 


0.157040 


100 


29.00826062334091 


29.16132759789657 


0.153067 


200 


58.16370616822418 


58.31479895068946 


0.151093 


300 


87.31958020522484 


87.47001658155840 


0.150436 


400 


116.4755621570468 


116.6256705858109 


0.150108 


500 


145.6315874033661 


145.7814989212340 


0.149912 



(b) 



TABLE I: We list ten typical values out of the whole database. 



The difference between AStoriji oo) and ASrecin ^ 
oo) is a function of m, AC(to) — Ctorirn) — Crecim). 
This difference is rooted in the topological difference be- 
tween a rectangle and a cylinder(since n — ^ oo, it is not 
a torus any more). The maximal value of AC(m) is 
AC7(2)=0.400162. AC(m) decays when m grows. The 
most rapid decay occurs from m — 2 to m = 50. This 
suggest that the topological difference between a torus 
and rectangle is most obvious in finite scale. When m 
continue to grow from m = 50 to m = 500, AC(rn) only 
drops from 0.157040 to 0.149912(Table|l|). 

The topological difference between rectangle and torus 
can be characterized by their different speed of entropy 
growth, 



FIG. 6: (a) The entropy growth of a rectangle from odd num- 
ber of columns to even number of columns, or vice verse. The 
Number of rows is fixed, m=10. (b)The entropy growth of a 
torus. The number of rows is kept to m=10. The number of 
columns varies from odd to even, or vice verse. 




AC(to) — AStor{n oo) — ASrec{n oo). (6) 

The rapid decay of AC(m) is illustrated in Fig. [51 Nu- 
merical check suggests that the decaying rate is much 
faster than exponential decay. The topological differ- 
ence between torus and rectangle will disappear when 
the width goes to infinity. However AC(m) will have 
finite non-zero value as long as m is finite. I call this 
property of AC(to) as topological stability. 
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FIG. 7: The limit value of entropy growth for n oo at 
different number of rows m. The square dot represents torus. 
The disc dot represents rectangle. The limit value of torus is 
always higher than that of rectangle. 



FIG. 8: The value of AC(m) = AStor{n ^ oo) - 
ASrec{n — > oo) as a function of the number of rows m. 
AC(m) approaches zero as m goes to infinity. 
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IV. ENTROPY REDUCTION BY CUTTING 
DIMER FILM ON FINITE RECTANGULAR 
LATTICE AND TORUS 

Experiment observation found two different liquid 
crystal molecules only dissolve in each other when they 
are separated by odd number of spacers [J]. If we con- 
sider the second law of thermodynamics, the fusion of 
two different types of molecule into one solution must 
increase the entropy of the system. The larger amount 
of entropy increasing, the easier for the two molecules to 
dissolve. According to this maximal entropy principle, 
one can calculate the entropy before the two molecules 
meet, Sa, and the entropy after they dissolve in each 
other, Sb- The value of entropy gain, AS* = Sb — Sa, 
determines wether the dissolve will occur or not. The 
practical liquid crystal system is too complicate to ac- 
complish an exact calculation of entropy. I only focus on 
the ideal theory of dimer model to show the entropy gain 
of fusing two rectangles with their length as odd number 
is much larger than fusing two rectangles with a length 
of even number. 

We take a rectangle of square lattice with m rows and 
n columns, and keep the number of rows m as one but di- 
vide number of columns n into two numbers: p and n—p. 
Then we calculate how much the number of dimer con- 
figurations increased when the two separated part fuse 
into one. 

As AZ is not always positive, if we take the logarithm 
function of AZ, negative AZ is excluded by the func- 
tion. Entropy is a monotonic function of the num- 
ber of states. When the number of dimer configura- 
tions increases (decrease), entropy will increase(decrease). 
Therefore we can use AZ to quantify the increase of en- 
tropy in the following. 

I made a sequence of cutting one rectangles into two 
daughter rectangles, {(1, n — 1), (2, n — 2), (3, n — 3), • • • }. 
The entropy gain for fusing two odd rectangles is much 
larger than that for fusing two even rectangles(Fig. [S] 
(a)). The entropy gain form two separated parabola 
band, the upper band is consist of all the divisions at 
odd numbers, the lower band includes all the divisions 
at even numbers. The same phenomena also occurs for 
cutting one torus into two daughter toruses(Fig. [9](b)). 
The entropy gain for fusing two odd torus is positive. But 
fusing two even torus does not increase entropy, instead 
the entropy after the fusion is reduced(Fig. [S] (b)). This 
is a significant difference from the case of rectangle lat- 
tice. The fusion of two rectangles obey the second law of 
thermodynamics, no matter their length is odd or even. 
The fusion of two odd rectangle is favored by entropy. 
But for torus, only fusing two odd torus will result in 
increasing entropy. 

The entropy gain for fusing two even toruses is negative 
at all different scales. I doubled the length of the mother 
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FIG. 9: (a) The increased number of states after dividing a 
finite rectangle into two smaller parts. The big rectangle has 
m = 18 rows and n = 18 columns, m = 18 is invariant. The 
length of n is divided into two parts, such as (1,17), (2,16), 
(3,15), • • • , and so forth. (b)The difference of total number 
states before and after dividing a finite torus with m — n — 36. 



torus by keeping m invariant. The gap of entropy gain 
between even case and odd case becomes smaller. The 
entropy gain of fusion at odd numbers approaches to zero 
in most case except the case that the cutting point is far 
from the middle. In the other case, I keep n invariant, but 
doubled m. The entropy gain for fusing two odd torus 
becomes zero. But entropy reduction for fusing two even 
torus becomes stronger. The gap of entropy gain between 
even and odd grows wider. Increasing n makes smaller 
entropy gap, while increasing m enlarger the entropy gap. 
Therefore the two spatial dimension parameter play two 
competing roles in controlling the entropy gap. 

In mind of the second law of thermodynamic, fusion 
of two rectangles prefer fusing two odd rectangles to fus- 
ing two even rectangles. While for torus, only fusing two 
odd toruses does not violate the second law of thermody- 
namics. In other words, a torus tend to divide into two 
even toruses instead of two odd toruses. Even though 
the over simplified dimer model is not directly related to 
the experiment observation of liquid crystal mixture, one 
can still get a rough primary understanding on why two 
polymers are miscible only when the linker includes odd 
number of units. 
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FIG. 10: (a) The variation of number of states for dividing a 
finite torus with (m — 72, n — 36) and uniting them. Com- 
paring with the finite torus in Fig. [9](b), the number of rows 
is doubled, (b) The difference of number of states before and 
after dividing the finite torus with (m = 36, n = 72). The 
number of columns is doubled comparing with Fig. [9](b). 

V. GEOMETRICAL AND TOPOLOGICAL 
CORRELATION OF TWO INTERSECTING 
LOOPS ON TORUS 




(a) (b) 



FIG. 11: (a) A torus with two topologically inequivalent cir- 
cles, Li and L2. (2) A torus covered by dimers can be con- 
structed from a rectangle lattice by taking periodic boundary 
condition in two spatial dimensions. 

From the point view of topology, a torus is a man- 
ifold generated by sweeping one circle around another 
circle(Fig. [TT] (a)). A torus is equivalent to a rectangle 
with periodic boundary condition. The two circles, Li 
and L2, are the boundaries of rectangle(Fig. [IT](b)). The 
entropy of dimer configuration is greatly reduced when 
the lattice manifold transforms from a torus into a rect- 



angle. This entropy reduction quantifies the geometric 
correlation of the two loops. 

In the free fermion theory of dimer model, the corre- 
lation of two monomers is defined by the ratio of total 
number of configurations with two monomer to that with- 
out two monomers[9] [10]. The two loops is a long string 
of many monomers. In most cases, the number of dimer 
configurations of rectangle is much smaller than that of 
torus. The subtle information of loop correlation lost in 
zeros. Thus I use a logarithm function to extract out the 
subtle information. The correlation of the two intersect- 
ing loops is defined as 

C(ii,i2) = log[(V'lV'2--"0nV'o'0lV'2---V'm)], 

= l0g[(iiL2)] (8) 

Loop Li is a string of ?i + 1 monomers, il'iip2 • • • tpni^o- 
Loop L2 covers m monomers, "i/;! t/;2 ••• V't?x, and the 
monomers at the intersecting point, V'o- This loop corre- 
lation function quantifies the correlation among (n -I- m -I- 
1) fermions. Following the monomer correlation in quan- 
tum field theory of dimer model, C(Li, L2) is essentially 
the entropy difference between rectangle and torus, 

L2) = log[|^] - Srec - Star. (9) 

When we keep total number of fermion as constant, 
the total length of the two loops is constant, Li -I- ^2 = 
n + m + 1 ~ const. The geometric correlation between 
the two loops can be understood by the shape of torus 
and rectangle. Shrinking Li will enlarge i2- If the circle 
Li of the torus becomes thinner, the other circle L2 must 
expands larger. For the rectangle, if Li becomes longer, 
L2 grows shorter. Thus geometric correlation between 
the two loops must be negative. 

The negative correlation of two loops shows different 
behavior according to even or odd number of fermion in 
each loop. First we keep m+n = 48 and both (rn,n) even. 
The correlation is weak at very small n and very large 
n, the maximal correlation is around the middle, n = 24. 
The opposite correlation behavior appears, when we still 
keep m + n = 48, but varies n and to in a sequence of 
odd number. The strong correlation now appears at very 
small n and very large n. Around n = 23 is the minimal 
correlation which is still larger than the maximal value 
for the case of even number sequences. One conclusion 
we can draw is the loop correlation for the case of odd 
number is much stronger than the case of even numbers. 
The correlation has a finite gap between odd and even 
case. The gap is closed when the length of the two loops 
grows to infinity. 

More over, we studied the correlation behavior for the 
first hybrid case that m is even and n is odd. The cor- 
relation behavior inherit half of the even case and half 
of the odd case. When m -I- n = 49, m is even and n 
is odd. The maximal correlation appear at small n, and 
the minimal correlation appears at large n. For the other 
hybrid case that m is odd and n is even, m -I- n = 49, 
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the correlation curve is the mirror image of the first case 
reflected by the axis n = 23. 



m 
n 



even 
■ even 



for L2 = 2A:. Then we define the correlation difference 
function as the difference between a correlation function 
with L2 = 2k + 1 and its neighboring correlation function 
with L2 = 2k, 

AC(Li, L2) = C(Li, L2 = 2fc + 1) - C{Li,L2 = 2k){10) 

If we calculate the limit value of AC{Li, L2) for the 
length of L2 = n growing to infinity, it leads to almost 
the same equation ^ in last section except a (—1) sign, 

lim AC(Li, L2) — ASrec{n 00) — AStorin oo](ll) 



(a) 



m — even 
n — odd 



Thus limLa^oo AC{Li, L2) = AC{m) is equivalent to the 
topological difference function of torus and rectangle de- 
fined in last section. Here AC(m) can be explained as 
the gradient of loop correlation. 



lim lim AC{Li,L2) = 0. 

Li— J-cxD L2— >-oo 



(12) 



This gradient function AC{Li,L2) has finite nonzero 
value as long as the length of the two loops arc finite. It 
is only when both the two loops become infinitely long, 
AC{Li, L2) will become zero. 



(b) 



m — odd 
n — odd 



VI. THE MELTING TEMPERATURE OF LONG 
POLYMER'S FILM 
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FIG. 12: (a)The sum of the two loops is constant, Li + L2 — 
49, i.e., m-|-n=48. Both m and n are even number varying 
from 2 to 48. (b) The total length of the two loop is 50, 
m-|-n=49. The length of L2 varies as following, (2,4,6,- ■ • ,48). 
The corresponding length of L2 is (47, 45, 43,- ■ • ,3,1). (c) The 
sum of the two loop is 49, m+n=48. Both m and n varies from 
odd number to odd number. 

Geometric correlation strongly depends on the scale of 
the two loops. While topological correlation is invariant 
when the two loops are expanding or shrinking. When 
the loop is expanding or shrinking, the total number of 
fermions along the loop will increase or decrease. Thus 
the correlation function for a fixed total number of par- 
ticles along the two loops is not topological invariant. 
However the difference between two selected correlation 
functions can be a topological invariant. We first fix the 
length of Li — m, and calculate the difference of two cor- 
relation functions one of which has odd number of length 
for L2 — 2A: -|- 1, the other has even number of length 



FIG. 13: Covering of a finite square lattice by long polymers. 
The length of the polymer showed above is 3, 4, 5, 6. The 
polymer represents the n-alkanes Cn-ff2n+2. Soft polymer can 
bend into different elementary configurations. 

The melting temperature for a film of long polymers 
can be computed by the same critical temperature Eq. 
([1]) for dimers. First we need to count the number of ex- 
ternal bonds. Then we count all possible ways to cover 
the lattice by tri-mer(3-connected monomers), quadra- 
mer(4-connected monomers), pentago-mer(5-connected 
monomers), and so on(Fig. [T3|) . However it is a difficult 
mathematical problem to exactly count all possible poly- 
mer coverings by long polymers. So far we only have the 
exact solution of dimer covering problem. If the polymer 
is rigid enough to keep the shape of a straight line, the 
entropy of longer polymer would be much smaller than 
the entropy of shorter polymers for covering the same 
area. We can always cut the longer polymers shorter, as 
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a result, the possible configurations would increase. Ac- 
cording to Eq. ([T]), smaller entropy induces larger melt- 
ing temperature. Although the number of external bonds 
also decreases for longer polymer, it falls far behind the 
decreasing speed of entropy. Quantitatively we can pre- 
dict that the melting temperature would increase when 
the length of the polymer increases. In mind of the ex- 
perimental observation of odd-even effect, we believe the 
exact entropy of long polymer covering would demon- 
strate odd-even effect. In fact, it is possible to measure 
the melting temperature of a film in lab. So experiment 
maybe is another way for solving mathematical problems. 

In reality, the n-alkanes C„i/2n+2 is not very rigid for 
large n. The soft chain would bend into different con- 
figurations(Fig. [T3l) . In that case, the entropy of longer 
polymer maybe is larger than shorter polymers. Then we 
would meet the anomalous odd-even effect. For a more 
complicate case, if we mix polymers of different length, it 
is almost impossible to get exact mathematical counting 
of all possible coverings. But we can divide the polymers 
into two classes: even-polymer and odd-polymers. The 
odd-even effect would still play a role. 

VII. SUMMARY 

Experiments observed that polymers with odd num- 
ber of length has better miscibility than polymers with 
even number of length. We believe mixing liquid crystal 
molecules of different length obeys the maximal entropy 
principle. Studying the dimer film can give us a basic un- 
derstanding on the odd-even effect. Since a long polymer 
of even number of length can be decomposed into many 
dimers. While an odd number of length is composed of 
many dimers and one monomer. 

We calculate the melting temperature and entropy 



growth rate of dimer film on a finite rectangle with fi- 
nite number of length. When the rectangle grows from 
an odd number of length to even number of length, the 
entropy increases much more than that for the case of a 
rectangle growing from an even number of length to an 
odd number of length. This results holds for all different 
initial length of the rectangle or torus. When the length 
of rectangle or torus grows to infinity, the speed of en- 
tropy increasing shows linear dependence on the width of 
rectangle or torus. The speed of entropy increasing for 
a torus is always higher than that of a rectangle. The 
difference between a torus and a rectangle decays to zero 
as the width grows to infinity. 

The correlation function of two loops on torus also de- 
pends on the odd-eveness of their length. We select two 
topologically inequivalent loops on torus. The two loops 
has one intersecting point and the total number of sites 
covered by the two loops is kept as constant. Then we 
enlarge or shrink one of the loops to odd number or even 
number of length. The correlation for odd numbers is 
much stronger than that for even numbers. 

We fuse two smaller rectangles into a bigger rectangle 
and compute how much entropy increased after the fu- 
sion. It suggest the entropy increasing for fusing two 
small rectangles with odd numbed of length is much 
larger than that for fusing two small rectangles with even 
number of length. While fusing two small toruses with 
even number of length into one big torus experience an 
entropy decreasing. The entropy increases only for fus- 
ing two small toruses with odd number of length. This 
odd-even effect is consistent with the second law of ther- 
modynamics. In a finite sized system, the entropy may 
experience small amplitude of decreasing periodically be- 
fore it finally reaches the maximal entropy. This entropy 
decreasing comes from the confinement of environment 
due to finite scales. 
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